Glassy behaviour in disordered systems with non-relaxational dynamics 
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We show that a family of disordered systems with non-relaxational dynamics may exhibit "glassy" 
behavior at nonzero temperature, although such a behavior appears to be ruled out by a face- 
value application of mean-field theory. Nevertheless, the roots of this behavior can be understood 
within mean-field theory itself, properly interpreted. Finite systems belonging to this family have a 
dynamical regime with a self-similar pattern of alternating periods of fast motion and trapping. 



The dynamics of disordered physical systems exhibits 
"glassy" features such as ergodicity breaking, slow dy- 
namics and aging. These systems are usually modelled 
by purely relaxational stochastic processes satisfying de- 
tailed balance. In this situation, glassiness arises from 
the complexity of the energy landscape: the representa- 
tive point strives to go downhill in energy along the maxi- 
mal slope while receiving random kicks from the thermal 
noise, and gets trapped into deeper and deeper energy 
valleys. 

Disordered physical systems have often been consid- 
ered a paradigm for complex behaviour in other fields, 
especially biology. However, there are no compelling rea- 
sons, in these fields, to restrict oneself to purely relax- 
ational dynamics. It is therefore important to know if the 
glassy properties exhibited by purely relaxational dynam- 
ical models are also present in the general case and, if so, 
whether the mechanisms responsible for glassy dynamics 
in non-relaxational systems are completely different 
from those acting in purely relaxational ones, even when 
the violation of detailed balance is small. If they were 
different, the picture that has evolved to explain the be- 
haviour of disordered systems would be entirely irrelevant 
as soon as non-dissipative forces are turned on. 

Hertz et al., Parisi || and especially Crisanti and Som- 
polinsky (CS) || showed some years ago that spin-glass 
behaviour was destroyed in several mean-field disordered 
models by an arbitrarily small but generic violation of 
detailed balance. 

We show in this paper that there is a wide class of 
disordered systems in which glassy behavior resists non- 
relaxational perturbations (even at T > 0), provided that 
either the initial condition is properly chosen, or that the 
size of the system is finite. (We shall not deal here with 
effects that are exclusive of the zero temperature case.) 
This behavior should be relevant for infinite systems in a 
finite dimensional space. The reason for this robustness 
can be understood, at least qualitatively, by a suitable 



interpretation of mean-field theory as follows. 

The stability of each energy (resp. free-energy) saddle 
can be characterized by the lowest eigenvalue A m i n of the 
Hessian. A class of models, of which that of Sherring- 
ton and Kirkpatrick (SK) is the best-known representa- 
tive, are marginal in the sense that in almost all states 
Amin goes to zero in the thermodynamical limit, and the 
single-state spin-glass susceptibility diverges. It is 

then not surprising that in such purely marginal mod- 
els a small (but still O(N)) perturbation may completely 
change their dynamics, as found in |^,^]. 

However, there is another class of models (currently 
thought to mimic 'fragile' structural glasses) having 
many non-marginal states with non-zero A m i n and for 
which the spin-glass susceptibility within such states is 
finite. Non-marginal states look locally much like fer- 
romagnetic states or the retrieval states in the Hopfield 
model. Now, it is easy to convince oneself that the sta- 
bility properties of a non-marginal state cannot be dra- 
matically altered by the combined effects of arbitrarily 
small non-relaxational forces and thermal noise fl(i[| . 

Consider for definiteness the typical case of the p-spin 
spherical spin-glass The energy landscape of the 

system fl3| features many saddle points of the energy 
function at different energy-density values £ < 0. One 
can identify a threshold value £th such that only the 
saddle points with £ < fth are minima. The lowest 
eigenvalue A m i n at each saddle point is proportional to 
the depth of the state beneath the threshold £th — £■ 
The only marginal saddle points are those just below the 
threshold, unlike the case of marginal models for which 
all states are marginal (or, in other words, for which there 
are no minima deep below the threshold) . From the rea- 
soning above, we might expect that only near-threshold 
minima will be destabilized by infinitesimal asymmetries, 
whereas the deeper a state, the more robust it will be. 

In the p-spin spherical glass, purely relaxational dy- 
namics starting from a random initial condition exhibits 
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"aging" |l4|: when N — oo, the system keeps touring 
a region just above the threshold, moving slower and 
slower but without ever getting completely trapped. For 
large but finite N (or for an infinite system in any finite 
dimensionality) the system penetrates a time-dependent 
amount below the threshold, and still ages due to the 
increasing depth of the traps it finds. In this purely re- 
laxational case, the dynamics is qualitatively similar to 
that of a marginal model. 

Now, just as a non-relaxational perturbation destroys 
aging in a marginal model, it also seems to destroy aging 
around the threshold in a non-marginal one. However, 
in non-marginal models, there are an infinite number of 
deeper states that remain stable in the presence of the 
perturbation. Hence, as soon as a finite system (or an in- 
finite system in finite dimensions) is able to penetrate be- 
low the threshold, it rediscovers the glassy features which 
had been destroyed above and near the threshold. 

In what follows we confirm this scenario for non- 
marginal models with non-relaxational dynamics. We 
first give evidence that mean-field dynamics (N = oo), 
starting from a random initial condition, yields for long 
times a time-translational invariant solution for the corre- 
lation and response functions even at small asymmetries 
(no aging), and that the correlation functions decay to 
zero (no ergodicity breaking) , confirming CS [^) . 

We then show, always within mean-field dynamics, 
that there are initial conditions such that the correlations 
do not decay to zero, in the presence of non-relaxational 
perturbations, even at T > 0. This confirms the exis- 
tence of stable regions with ergodicity breaking. Obvi- 
ously, such regions are never found if the system is infinite 
and starts from a random configuration. 

Finally, in order to estimate the time scales involved 
for trapping in a large but finite system, we performed 
simulations. 

We consider a system of N variables s = (s±, . . . , s;v), 
subject to forces F t given by 

F l (s)= ■tf W *" 1 «*-*J,- 1 . (!) 

{ii,—.3*>-i} 

where the couplings are random Gaussian variables. For 
different sets of indices {i, ji, . . . , j p -i} the J's are uncor- 
rected, while for permutations of the same set of indices 
they are correlated so that 

Fi{s')Fj(s) = SMq) + s' iSj f 2 (q)/N, (2) 

where q = (s • s')/N. In the purely relaxational model 
one has f 2 (q) = f[(q)- We consider here f 2 (q) = af[(q), 
where fi(q) = pq p_1 /2. The purely relaxational case has 
symmetric couplings j^ 1 '" Jp ^ 1 under i <-> jk (a = 1) while 
for uncorrected jf - jp - 1 and jf 2 - j "- 1 : a = 0. 

We consider: (i) the continuous spherical model s| 2 = 
N, with Langevin dynamics s\ = — F,;(s) — z(t)si + 



hi(t) + T)i(t), where rj is a white noise of variance IT — 
2//3, z(t) is a Lagrange multiplier enforcing the con- 
straint and hi(t) is an external field (usually set to 
zero); (ii) the Ising s$ = ±1 model with Metropo- 
lis dynamics, in which a randomly chosen spin flips 
with probability min[l, exp(— 2f3siFi)]. The "energy" 
£ is defined in all cases by the expression N£(s) — 
— (F • s)/p to which only the symmetric part of the 
force contributes. We monitor the autocorrelation func- 
tion C(t,t') = (s(t) ■ s(t'))/N and the response function 
G{t,t>) = Z l 5{s t (t))/Sh l (t')/N. 

Both models are marginal for p — 2. The dynamics 
of the asymmetric p = 2 spherical and ±1 (SK) models 
were studied numerically and by mean-field theory (large 
N) by CS, who found no glassiness at finite T fl5| . Here 
we concentrate on the non-marginal case p > 2. 
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FIG. 1. Interruption of aging: C(t,t') vs. t — t' for 
t = 200,400,600,800 at T — 0.2 and a = 0.96 compared 
to aging for a = 1 (inset). 

We have solved numerically the mean-field equations 
(exact for N — > oo) for the correlation and the response 
function of the spherical (p = 3) model, starting from a 
random initial configuration [^6|. The result for a = 0.96 
and T = 0.2 is shown in Fig. 1. The corresponding curves 
for a — 1 are shown in the inset. On can see that there is 
an initial regime in which the solution exhibits aging as 
for a = 1, followed by a crossover to time-translational 
invariance (no aging). The closer a is to one, the longer 
the (interrupted) aging regime lasts. We have not found 
any evidence for a critical asymmetry (a c < 1) beyond 
which the system ages forever although we cannot rule 
out a transition for a even closer to one. The correlations 
decay to zero ruling out ergodicity-breaking also as in CS. 

In order to prove that there are trapping regions for 
a < 1, T > 0, we have studied the mean- field dynamics 
starting from a configuration with given (low) energy, 
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but otherwise random jL2]. This amounts to solving a 
static problem for the initial condition, at a temperature 
T' tuned to give the desired value of the energy Jl7| . 
Figure 2 shows that for a > 0.86 and sufficiently small 
initial energies the system remains trapped at non-zero 
temperature (forever, if JV = oo). 
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FIG. 2. C7(i w + T,t w ) vs. r for t w = 100 at T = 0.1 
starting from an initial condition with low energy-density 
(X" = 0.1). From top to bottom the asymmetry parameter 
equals a = 0.88, 0.86, 0.84, 0.82, 0.8. 
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FIG. 3. The magnetization m(t,t') = J* dsG(t,s) vs. 
C(t, t') starting from a random initial condition for T = 0.2 
and q = 0.3, 0.7, 0.9. Inset: same, starting from low energy, 
T = 0.1, a = 0.8,0.9. 

Figure 3 shows the fluctuation-dissipation (FD) ratio 
G(M') = X[C(t,t')](9(t-t')/T)dC(t,t')/dt' Q starting 
from random initial conditions. The same plot for a — 1 
consists of two straight lines corresponding to X = 1 for 



C > qsA (FDT) and X = constant < 1 for C < qsA- 
It is remarkable that the plot crosses over smoothly, as 
a — ► 1, to the one holding for the relaxational case (note 
that X[C] < 1, VC, Va). In the inset we plot the FD 
ratio starting from low energy, for a = 0.8 (untrapped) 
and a — 0.9 (trapped). 

In conclusion, the mean-field analysis reveals that 
the present model has trapping regions that cannot be 
reached when N = oo because the time needed to fall 
into them diverges with N and, by the same token, that 
the time needed to escape from a trap also diverges. 
The question how typical falling and escaping times scale 
with N is beyond the present analytical tools. We have 
thus performed numerical simulations, choosing for con- 
venience the ±1 version of the (p = 3) model. 

Figure 4 shows the typical behaviour at T — 0.01 and 
a = 0.5, for N = 50. The system alternates between 
periods of trapping and periods of rapid motion at high 
energy. A blow-up in time of the same run shows that 
the overall appearance of the graph is self-similar. The 
longest trapping time is of the order of the total observa- 
tion time, which indicates a broad distribution of release 
times. We have found such behaviour, for N = 50, in a 
region in the (T, a) plane bounded by a ~ 0.4 at T ~ 
and T = T c ~ 0.05 at a = 1. Traps are visited once show- 
ing that there is a large time-span between the smallest 
falling time and the maximal ('equilibration') trapping 
time. For larger system sizes, N = 100, 200 falling and 
escaping times increase with N, as expected. 

This behaviour is reminiscent of the non-relaxational 
dynamics of a particle in a random velocity field ]22] ] 
where anomalous diffusion is due ]2^,^3j to broad, Levy- 
stable distributions of trapping times. Indeed, this model 
likely to be close to a microscopical realisation of the re- 
lated 'trap model' that has been fruitfully used to de- 
scribe aging in [Q . 

We have therefore shown that aging and ergodicity 
breaking resist non-relaxational perturbations in the dy- 
namics of disordered systems if they have non-marginal 
states. Finite-size fully connected systems (relevant, e.g., 
for modelling biological networks) may thus exhibit strik- 
ing aging effects. This is also likely to be the case for 
finite-dimensional systems whose mean-field limit is non- 
marginal. 

Let us also mention that we have checked that sys- 
tems with random forces deriving from a potential and 
strongly perturbed by random non-potential forces cor- 
related with a different range exhibit aging phenomena 
even at the mean-field level. 

Further scenarios can be envisaged: it is likely that 
some no-go results obtained almost a decade ago did not 
exhaust all possibilities of nature and that this direction 
is open for further research. 
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FIG. 4. a. Instantaneous energy £ vs. time (Monte-Carlo 
sweeps) in the ±1 (p = 3) model for a — 0.5 and T = 0.01. 
b. A blow-up of the first 10 000 sweeps. 
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